Field dependent dynamics in the metallic regime of the half-filled Hubbard model. 
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A systematic study of the effect of magnetic field (h) on Hubbard model has been carried out 
at half filling within dynamical mean field theory. In agreement with previous studies, we find a 
zero temperature itinerant metamagnetic transition, reflected in the discontinuous changes in mag- 
netization as well as in the hysteresis, from a paramagnetic (PM) metallic state to a polarized 
quasi-ferromagnetic (QFM) state, at intermediate and large interaction strength (U). The jump 
in magnetization vanishes smoothly with decreasing interaction strength, and at a critical U, the 
transition becomes continuous. The region of 'coexistence' of the PM and QFM solutions in the 
field-C/ plane obtained in this study agrees quantitatively with recent numerical renormalization 
group calculations, thus providing an important benchmark. We highlight the changes in dynamics 
and quasiparticle weight across this transition. The effective mass increases sharply as the transi- 
tion is approached, exhibiting a cusp-like singularity at the critical field, and decreases with field 
monotonically beyond the transition. We conjecture that the first order metamagnetic transition 
is a result of the competition between Kondo screening, that tries to quench the local moments, 
and Zeeman coupling, which induces polarization and hence promotes local moment formation. A 
comparison of our theoretical results with experiments on ''He indicate that, a theory of '^He based 
on the half-filled Hubbard model places it in a regime of intermediate interaction strength. 



PACS numbers: 

I. INTRODUCTION 

One of the most interesting aspects of strongly cor- 
related electrons is the self-consistent emergence of a 
low energy scaled. This has been demonstrated over the 
years from a large number of theoretical and experimen- 
tal studies. The Kondo quenching of local moments and 
the emergence of a Fermi liquid below a critical tem- 
perature in Mott-Hubbard systems is clearly seen in dy- 
namical mean-field theory (DMFT). A question, there- 
fore, naturally arises as to how an external magnetic field 
would interact with the local moments and how the Mott- 
Hubbard physics gets affected by this added field that 
couples with the moments. 

The Mott-Hubbard systems Hkc {Vi-xCr^)20^, 
V20^ display local magnetic moments which also inter- 
act through residual antiferromagnetic exchange. There 
is a competition between the alignment of the local mo- 
ments in presence of a magnetic field and the residual 
antiferromagnetic exchange. Magnetic transitions are 
known to occur in the half filled Hubbard model^t'^i in 
the insulating state. Of course, the problem is more 
involved in the metallic side due to the presence of an 
additional energy scale (the local Kondo temperature), 
associated with the quenching of local spin fluctuations 
at low temperature. So, the question raised above is 
much more subtle in the metallic phase of the Mott- 
Hubbard systems. A quantitative study of the dynam- 
ics and transport properties of this model is not easy. 
Recently, much progress has been made with the ad- 
vent of dynamical mean field theorj*^^— which is exact 



in the limit of infinite spatial dimension. Within DMFT, 
a correlated model maps onto a single impurity model 
in a self-consistent conduction-electron bath. Meth- 
ods like Numerical renormalization group (NRG)^, Lo- 
cal moment approach (LMA)ii"— , Exact diagonalization 
(ED)I, quantum Monte Carlo (QMC)^, Iterative pertur- 
bation theory (IPT)ii"— have been used to study the 
impurity problem within the DMFT framework. Mott- 
transition has been observed in the half filled Hubbard 
model using DMFTii^. 

Prior to the advent of DMFT, there have been two 
major approaches for the strongly correlated Fermi liq- 
uid - the Stoner approacbi^"— and the Gutzwiller ap- 
proach (GA)^"— . The former approach views the sys- 
tem as close to a ferromagnetic transition and the lat- 
ter views it as a strongly correlated, nearly incompress- 
ible Fermi liquid close to a Mott localization. While 
the Stoner theory is essentially a high temperature ap- 
proach, the Gutzwiller method starts from a Hubbard 
model description of the strongly correlated electrons 
and projects out the local double occupancies from the 
ground state in a variational calculation. The GA is in- 
capable of describing the metallic phase of the half-filled 
Hubbard model, and also fails to incorporate any dy- 
namics or temperature dependence. However, it docs 
have the merit of being analytically tractable, and thus 
can yield qualitative insight, provided the results have 
been benchmarked against other detailed calculations. 
These two approaches have very different predictions for 
the response to an external magnetic field. The Stoner 
approach predicts a smooth variation of magnetization 
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in presence of a magnetic field while the Gutzwiller ap- 
proach gives a first order transition of the magnetization. 
Experimentally, an effort to distinguish between the two 
predictions, based on studies in liquid ^He^, finds a 
smooth variation of the magnetization as predicted in 
Stoner approach. Within DMFT, calculations^^ using 
exact diagonalization (ED) as the impurity solver were 
carried out to study the effect of magnetic field on the 
Hubbard model. A fictitious temperature scale was intro- 
duced through a frequency cutoff. The calculations were 
carried out at a fixed magnetization. A metamagnetic 
transition was obtained, and the T — U 'phase diagram' 
was extended to the finite field axis, that demonstrated a 
similarity between the effects of temperature and field. A 
qualitative agreement with the Guzwiller approach was 
also found. Later, an exchange term was introducedS^ to 
incorporate the spin-spin correlation to combine features 
of Mott localization and almost ferromagnetic Stoner in- 
stability. Changes in dynamics with applied field, how- 
ever, was not studied in these papers. A recent study of 
the same problem was carried out by Bauer— using the 
numerical renormalization group method within DMFT. 
The focus in this study was on the metamagnetic tran- 
sition, and the results were qualitatively similar to the 
earlier work; the main result being that the divergence 
of susceptibility at the metamagnetic transition is not 
a consequence of effective mass divergence, but occurs 
through the quasiparticle interaction term. The predic- 
tion for the coexistence region in the h — U plane was not 
made explicitly in the NRG study, nevertheless, it is pos- 
sible to deduce this, from the results presented. There 
appears to be a strong disagreement between the NRG 
and ED study in this context. 

In this paper, we study the effects of magnetic field on 
correlated electron liquids within the framework of the 
half-filled Hubbard model, using the iterated perturba- 
tion theory (IPT) as the impurity solver for the DMFT. 
The relative merits and demerits of the IPT are discussed 
in the next section. Since DMFT is known to treat the 
local dynamics quite accurately, it is expected to reveal 
the local nature of the competition between the spin fluc- 
tuation and the aligning field. We calculate the single 
particle dynamics, magnetization and other properties of 
the Hubbard model in the presence of an external mag- 
netic field. In the infinite dimensional limit, the mag- 
netic field appears only as a Zeeman term and there is 
no explicit orbital contribution in the Hamiltonian as the 
lattice model is mapped onto a single-impurity model. 
There are no non-local terms, and the effects due to 
the non-interacting medium are included as a dynamical 
mean field in the hybridization. Our focus in this study is 
on the field-induced changes in quasiparticle weight and 
spectral functions as the interaction strength is varied. 
Although the DMFT-hNRG study by Bauer^Si does ex- 
amine dynamics, the detailed behaviour of the crossover 
from the low field Fermi liquid behaviour to the high-field 
non-Fermi liquid behaviour was not examined. Such a 
crossover was proposed in the DMFT+ED work^^, but 



was not supported by dynamics. Here we provide a com- 
prehensive picture of this crossover. In our study, we 
show the detailed behaviour of the quasiparticle weights, 
that show an asymmetric cusp like singularity in the ef- 
fective mass as a function of field. We show that the 
effective mass increases monotonically as the metamag- 
netic transition is approached, and beyond the transition 
the effective mass decreases smoothly and monotonically 
for lower U, while for higher U, there is a sharp discon- 
tinuous decrease. The applied field gets renormalized to 
an effective field through the polarization of the medium 
(within DMFT) and via interaction effects. 

We show that the lattice Kondo resonance at the Fermi 
level splits into two at large fields: the up and down com- 
ponents move away from the Fermi level and finally form 
a spin polarized band insulator. The shift of the band 
is not rigid due to the competition between the Kondo 
screening and the Zeeman effect. The hysteresis across 
the metamagnetic transition is used to find the 'phase 
diagram'. Our calculations agree quantitatively with the 
NRG study of Bauer— and we regard this agreement as 
an important benchmark. The low field universal fea- 
tures of the spectral functions are carefully examined 
and adiabatic continuity to the non-interacting limit is 
demonstrated. We argue that the field-induced meta- 
magnetic transition is from a correlated paramagnetic 
metallic liquid to a band insulator, and in this sense, is 
very different from the interaction driven first order Mott 
transition, which is from a correlated metallic liquid to 
a correlated Mott insulator. We compare our theoretical 
result with the experimental one on liquid ^He and find 
good qualitative agreement. 

The paper is organized as follows. We begin in Sec. II 
with a brief description of the model. In Sec. HI, we 
present our theoretical results and their analysis. We 
also compare our theory with the experimental results 
on liquid ^He. Finally, we conclude in Sec. IV. 

II. MODEL AND FORMALISM 

The single band Hubbard model Hamiltonian describ- 
ing correlated electrons in the presence of an external 
magnetic field (the orbital contribution is neglected) is 
given by: 

■1(7 (1) 

The first term describes the kinetic energy of the non- 
interacting conduction (c) band due to nearest neigh- 
bor hopping t. The second term refers to the on-site 
repulsion U . The final term represents the orbital en- 
ergy and the Zeeman splitting in an external magnetic 
field. For the particle-hole symmetric case (half-filled) 
case considered in this work, the orbital energy is given 
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by ec = —U/2. In the limit of large dimensions, D — )■ oo, 
the hopping needs to be scaled as t oc t^,/^/15. We choose, 
for convenience, a hypercubic lattice, for which the non- 
interacting density of states is an unbounded Gaussian 
(Po(e) = exp(-e^/tj)/(^t*)). We set = 1 in the 
following calculation. 

Within dynamical mean field theory (DMFT), which 
is exact in the limit of infinite dimensions, as is well 
known, the lattice model may be mapped to an effective 
single-impurity Anderson model with a self-consistent hy- 
bridization. The major simplification that occurs is that 
the self-energy and the vertex function become local and 
momentum-independent. In the presence of an external 
global magnetic field, the local, retarded Green's function 
is given by. 



G„{u,h) = HYl„{uj-h)] 



(2) 



where 7ct(w; h) = w+ — Cc + fft, — Yj„{uj, h), H[z] is the 
Hilbert trasform defined as 



H[z] 



and Scr(w; h) is the spin-dependent local self-energy. The 
local Green's function may be used to define a hybridiza- 
tion function S{lo] h) as 



G.(w,/i) 



1 



7cr(w; h) ~ S{u}; h) 
or, expressed another way, 

1 



(3) 



(4) 



The hybridization function defines the host or the 
medium within which the impurity is embedded, and the 
non-interacting host Green's function is then given by 



Gaiuj,h) 



1 



(5) 



a;+ — Ec + crh — S{uj: h) 
which is equivalent to using the Dyson's equation 

g~^{uj,h)^G-\uj,h) + j:^{uj,h). (6) 

This host Green's function is then used to construct a 
new impurity self energy, which if substituted in equa- 
tion [21 yields a new lattice Green's function. Thus the 
solution within DMFT proceeds most easily in an itera- 
tive manner using equations (2)-(6) until self-consistency 
is achieved. 

Although easily stated, constructing a new impurity 
self-energy for an arbitrary non-interacting host has re- 
mained a major bottleneck in the DMFT scheme. Var- 
ious impurity solvers have been developed, most impor- 
tant of which are Quantum Monte Garlo, exact diag- 
onalization, numerical renormalization group, diagram- 
matic perturbation theory based approaches etc. We em- 
ploy the iterative perturbation theory (IPT) approach to 



solve the impurity problem within dynamical mean field 
theory (DMFT). The IPT is a very simple yet power- 
ful diagrammatic perturbation theory based approach. 
It has been benchmarked extensively against more ex- 
act methods in the zero field half-filled Hubbard model 
problem. The agreement, though qualitative, is excel- 
lent. The main advantage of IPT is the ability it provides 
to capture real-frequency dynamics, and hence trans- 
port on all frequency, temperature scales and interaction 
strengths. For a more detailed comparison with other 
impurity solvers, we refer the reader to a previous article 
by one of the authors of the present worki^. 

The zero temperature (T = 0) self energy has the fol- 
lowing form, consisting of the static Hartree contribution 
and the dynamical part. 



(7) 



where n„ = I?cr(w) and ^^{uj) = —^lmQa-{Lj) is 

the spectral function of the host Green's function. The 
local IPT self energy (w, h) which satisfies the Lut- 
tinger theorem automatically at half-filling is given by 
(this zero field algorithm has been discussed recently in 
detaiiil). 



E(,2) {uj,h) = U^ duj'Vg (-w; h) X 

J —OQ 

[xfioJ + oj')f{u:') + xr(w + ^')/(-^')] (8) 



where 



xfH= / dL0'V,{J)Vg{L0+J)f{-L0')f{uj+L0') 

J —oo 

(9) 

The magnetization can be calculated from 

m = (nt) - (n;) (10) 
where (?v) = /"^ ^)lmG „ {uj) . 

III. RESULTS AND DISCUSSIONS 

In this section, we will discuss the results of our calcu- 
lations so far on the effect of magnetic field on Hubbard 
model at half filling. Before discussing the finite field 
results, we will briefly review the well-known zero field 
results for the Hubbard model. At zero field, the system 
goes from a metal to an insulator at U > Uc^- Below 
Uci{< Uc2)i the system is metal and shows Fermi liquid 
behaviour. There is a coexisting region between Uci and 
Uc^ . The values of C/ci and Uc^ for the hypercubic lattice 
are 3.67t* and 4,.l%t:^^. We are interested in the region 
below [/c2- At zero field, the spectra shows an universal 
scaling form in terms of (where lol = Zi* and Z is 
the quasiparticle weight). 
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We will begin with a discussion of the symmetry prop- 
erties of the Green's function and the associated self en- 
ergies. The following symmetry holds in the half- filled 
case: 

G„{w;h)^-[G^{-iu;h)]* (11) 
Ea{oj;h) = -[j:si^uj;h)]* (12) 

To see this, imagine carrying out the self-consistent 
DMFT iterations beginning with a symmetric hybridiza- 
tion function, as appropriate for the zero-field case. 
Naturally, the resulting host Green's function (equa- 
tion [5])will satisfy the above symmetries. The self- 
energies calculated through the IPT ansatz (equations [8]- 
in]) would then also satisfy the same symmetries, as then 
would the full interacting Green's function (equation [3]) . 
One consequence of such a symmetry is that the quasi- 
particle weights Z^^ would be spin- independent. 

Rewriting equation [7] as T,a{<^',h) = ^{h — am) + 

Eo- (w; h) where fh = hf — is a, fictitious local mo- 
ment (of the host/medium); and using ec ~ —U/2 and 
n = 1, we transform Eq. [2]as, 



J-oo UJ+ — e + cr (/i + ^nij — Y!i^\uj\ h) 

(13) 

We would like to mention a technical point here. We 
found that the IPT equations arc ill-behaved for a fixed 
field. In other words, if we wish to find the value of 
the magnetization self-consistently through the iterative 
procedure, keeping the field fixed, we run into problems 
of convergence. An easy way to avoid this instability is 
to work with a fixed value of the composite number hm = 
h + Ufh/2. The equations converge easily and a unique 
solution is found. The host Green's functions can be then 
used to find to, and the field value is obtained a.s h ~ 
hm — Uffi/2. A similar strategy was employed by Laloux 
et al2^ in their ED+DMFT approach, where they had to 
use a fixed magnetization to achieve convergence. The 
final value of the actual magnetization is found through 
equation [TUl 

Our results for magnetization are summarized in fig- 
ure [TJ The magnetization (to) is shown as a func- 
tion of field {h) for increasing values of the interac- 
tion strength U /<* . The magnetization smoothly ap- 
proaches the saturation value of 1 for J7 < 1.4t*, while 
for U > 1.4i*, a metamagnetic first order transition is 
seen. The magnetization jumps at a finite value of the 
field hci which falls sharply with increasing [7/^* (fig- 
ure |3]), and beyond C//i» ~ 2.2 becomes unresolvably 
small within the present approach. The jump size in- 
creases from at a Uhc = 1.4t* to a large value (^ 0.67) 
at U = 2. It*. The value of Uhc = 1.4^* does agree 
quantitatively with that obtained by Bauer— through 
NRG+DMFT. In the latter, a semi-circular density of 
states has been chosen to represent the bare conduction 
band, pI""' = 2VD^ - e^/nD^ with D = 2U. The value 




Field (h) 

FIG. 1: (Color online) Magnetization as a function of field for 
various U/t, values. The magnetization smoothly approaches 
the saturation value of 1 for U < lAt,, while for U > lAt,, 
a metamagnetic first order transition is seen. The magneti- 
zation jumps at a finite value of the field, which falls sharply 
with increasing U/t* (see figure (31), and beyond U/t, = 2.2 
becomes unresolvably small within the present approach. The 
inset shows the same data as the main panel on a logarithmic 
field axis, so that the sharp decrease in the he value is seen 
clearly. 



of Uhc/t* was found to be 2.64 in the NRG study22, which 
corresponds to a Uhc/D ~ 1.32, and thus the agreement 
with the IPT value (1.4) for the same provides our results 
an important benchmark. (1.4) This metamagnetic tran- 
sition has been observed and studied previously within 
the GA and DMFT approaches. The low field state is 
a paramagnetic state, while the high field state beyond 
the transition is a polarized quasi-ferromagnetic state. 
As with the zero field case, the metamagnetic transition 
being first order, would be accompanied by a hystere- 
sis, as the field is swept through a closed cycle across 
the transition. This is shown in figure [2] for U = 1.5i,. 
Thus, there exists, for a range of U values, a region in the 
hc — U plane that defines the 'coexistence' region, i.e the 
region, where the paramagnetic metallic and the quasi- 
ferromagnetic insulating solutions coexist. This region 
may be identified through the bounds of the hysteresis 
curve for each U. In figure [31 the shaded region repre- 
sents the coexistence region, and the solid lines are the 
equivalents of Ud and Uc2 for the finite field case. In- 
deed, Bauer— reports that the metallic solutions were 
impossible to get for U/t^, > 4.5, which corresponds to 
U / D c:i 2.2, which is the right boundary of the coexis- 
tence region in figure [Sj and thus the coexistence region 
obtained within IPT for the hypcrcubic lattice agrees ex- 
cellently with that of the NRG study. Since the quasi 
ferromagnetic high field insulator is not however a Mott 
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FIG. 2: Hysteresis is seen across the metamagnetic transition. 
We illustrate this for U = 1.5t* here. The dotted line is the 
actual obtained line when the parameter hm is continuously 
increased. The region where the susceptibility, x = dm/dh 
becomes negative is excluded to obtain the hysteresis. 

insulator, the term 'equivalent' above is not strictly pre- 
cise. In fact, this two-phase region at finite field is not 
an "extension" of the Uci — Uc2 region in the ft, = 
plane onto the finite h plane. They do not merge onto 
each other and represent very different spinodal regions 
across different states. 




FIG. 3: The coexistence region (shaded) as obtained within 
IPT for the hypercubic lattice. At large U{> 2.2t,), the sys- 
tem is susceptible to being polarized for an infinitesimal field. 

We shift our focus to the quasiparticle weights or in- 
verse effective mass. In figurelU we show the Z as a func- 
tion of field for various U values, ranging from low U at 
which metamagnetism is weak or absent, to a relatively 
high U, i.e the strongly metamagnetic region, where the 
system is susceptible to being polarized even for an in- 
finitesimal field. 

At low U = 1.0, the magnetization increases smoothly 
with increasing field, and the effective mass decreases 
monotonically, approaching the bare mass at high fields. 
The behaviour of m and Z changes qualitatively, for 
U ~ 1.3t*; the effective mass displays a cusp, while 
the magnetization changes slope significantly. The cusp 
sharpens for U = 1.5^*, while Z suffers a discontinuous 
increase at the metamagnetic transition. At even higher 
U = 2.0, both the magnetization and effective mass dis- 
play discontinuous changes at the transition. Thus the 
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FIG. 4: The dependence of quasiparticle weight (inverse ef- 
fective mass) is illustrated here for four U values. The cusp 
singularity is seen to develop for intermediate and high U, 
and at precisely the metamagnetic transition. 

field-dependence of magnetization and effective mass is 
very sensitive to the specific U/t^ values under consider- 
ation. 

The main advantage with the IPT approach used here, 
is that we can study the real frequency dynamics in de- 
tail. We now proceed to elucidate the field dependence of 
the dynamics before and across the metamagnetic transi- 
tion. We begin by exploring the universal scaling proper- 
ties by carrying out a low frequency Fermi liquid analysis. 
In the metallic phase of interest, the real part of self en- 
ergy may be expanded about the Fermi level to first order 
in oj, 

Re42)(^,/i) ^ Re42)(o,/,) + (1 _ (14) 

Za[h) 

where Z„{h) = [1 - dJ:^{uj; h)/duj]-^. For the half-filled 
case considered in this paper, Z-^ = Z^ (as argued before), 
hence we drop the spin subscript. The imaginary part of 
the self energy remains Fermi liquid like (^ 0{uj^)) be- 
cause of the structure of the IPT equations (equations [S]- 
[5]). A finite imaginary part can be acquired by the self 
energy only at finite temperatures in this approach. Sub- 
stituting the above in equation 1131 the spectral function 
D„ — — ImGcr/TT, is obtained as 

D^{u:,h)'^p^{Cj + ah'ff) (15) 

where w = lo/wl (with wl ~ Z{h)t^) and the spin- 
independent (because of the symmetries, equations I12[) 
effective field is given by 

h^ff = h+ j7fi - (TReS(,2) (0, h). (16) 

The above Eq. [15] is the renormalized non-interacting 
limit of the full spectral function. Given the above form 
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FIG. 5: (Color online) The effective field is plotted against 
applied field (h) for the U values chosen in figure [l] 

of the spectral function, it is easy to see that adiabatic 
continuity to the non-interacting limit is achieved pro- 
vided the real part of the self energy may be expanded as 
in equation 1141 Using eauation llSl we can infer that the 
spin-summed spectra would also be a universal function 
of w for a fixed /le//- Whether the full field-dependent 
spectra satisfy adiabatic continuity and scaling may be 
tested through equation [15] 

The dependence of effective field (ft-e//) on the applied 
field is shown in figure [5] for various U values. The 
metamagnetic transition is reflected in the effective field 
as well. A careful examination of the contribution to the 
effective field indicates that the discontinuity is present 
not only in the real part of S'^' (0; h), but also in the host 
magnetization fh. 

The spectra for various U /i, plotted clS cL function of 
lo/ojl = u!/{Z{h)tf,) collapse in the neighbourhood of the 
Fermi level as shown in figure |6l Also shown in the same 
graph is the non-interacting density of states (dos) for 
a field of /i = 0.3i, and Z = 1. All the [/ > spec- 
tra are seen to be identical to that of the non-interacting 
case, near the Fermi level, thus exhibiting adiabatic con- 
tinuity and universal scaling as a function of lo/ujl, and 
validating equation 1151 This demonstrates Fermi liquid 
behaviour for all h lower than the for a given U /t* . 

We now move on to the field dependence of the spectral 
functions for a given interaction strength. In the absence 
of clear metamagnetism ([/ = 1.0, the spectral function 
evolution with field is very similar to a non-interacting 
case, as expected. The density of states at the Fermi 
level continues to decrease with increasing field, until the 
up and down spin bands move so far from the Fermi 
level, that a band insulator is obtained. This is shown 
in the upper panel of figure [T] The lower panel of the 
same figure represents U = 1.3t*, which is close to the 
onset of metamagnetism. Again, the behaviour is sim- 
ilar, albeit the decrease of the density of states at the 
Fermi level with increasing field is more pronounced. We 
add here that the IPT results for the field-dependence 
of the spectral functions are not in agreement with the 



FIG. 6: (Color online) Adiabatic continuity and scaling: The 
full spectral function when plotted as a function of (jj/ojl = 
uj / {Z(h)tt) for various U (shown next to the curves) collapse 
in the neighbourhood of the Fermi level. The universal form 
is seen to be identical to the non-interacting spectral function 
for the same parameters (see text for discussion). 
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FIG. 7: (Color online) The upper (lower) panel represents 
the spectral function for U — 1.0 {U — 1.3) and various field 
values plotted against the bare frequency, uj/t,. 



NRG-I-DMFT resuhs of Bauer et. al2i. In the latter, the 
spectral function seems to be pinned at the Fermi level 
until the metamagnetic transition, whence there is a sud- 
den drop in the density of states. We do not understand 
this difference with NRG completely, but speculate that 
this disagreement could be due to the non-conserving na- 
ture of IPT. 

The evolution of the spectral function for [/ = 1.5 and 
U = 2.0 (figure [5]) with field is quite dramatic, reflect- 
ing the metamagnetic transition. The spin bands shift 
rapidly with field, and for a critical field at which the 
transition happens, the shift from a metallic to a quasi- 
ferromagnetic insulator takes place. Though the insulat- 
ing state is reminiscent of a Mott insulator, it is actually 
far from being one. The two 'Hubbard bands' are just 
the two spin bands shifted away from the Fermi level, 
and broadened to an extent that they just resemble non- 
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FIG. 8: (Color online) The upper (lower) panel represents 
the spectral function for U = 1.5 (U = 2.0) and various field 
values plotted against the bare frequency, oj/t*. 



interacting Gaussian density of states. In fact, the large- 
field spectral function may be deduced from the large 
field asymptote of equation [131 At large field /i ~ , 
the effective field becomes comparable to or larger than 
the largest scale in the problem, i.e U in strong cou- 
pling. This is because, for h t^,, the magnetization 
saturates, m — > 1, hence the effective field becomes at 
least h + U/2. The real part of the second-order self 
energy contributes 1/U aX frequencies a; ~ f/, because 
S(^)(u;) ^ l/w at large frequencies. Thus from equa- 
tion [13] the large field spectra must be simply given by 
[Pq{ijJ + heff) + pq(u} — h(.ff)] /2. This is just a sum of 
two Gaussians each centred around ih^ff respectively. 
Indeed, from figure [5l we see that for U = 2.0, the value 
of the effective field is ^ 2.5 for a field of /i = 0.5; and 
the lower panel of figure [5| shows that the two 'Hubbard 
bands' are centred around ~ ±2.5. This indicates, that 
at high field, all correlation effects are lost, and the 'Hub- 
bard bands' are nothing but the bare density of states 
shifted away from the Fermi level, and thus the insulator 
at high fields is just a band insulator. 

The split between the up and down spin bands may 
be quantified, and is shown in figure [Sj for various U val- 
ues. The field induced shift of the spin bands is seen to 
be highly non-linear and indeed even discontinuous for 
U ~ 1.5 and U = 2.0, naturally at the metamagnetic 
transition. In the non-interacting case this split is a rigid 
Zeeman shift. But in the interacting case the shift is not 
rigid at all. In the periodic Anderson model (appropriate 
for heavy fermion systems), the effect of magnetic field 
in the Kondo lattice regime is similar (without metamag- 
netism), and it has been argued that the non-rigid shift 
is due to the competition between Kondo screening of 
the local moments (whose effect is to quench the local 
moments) and the Zeeman coupling (whose effect is to 
polarize the local moments). In the present case, since 
the lattice Hubbard model is mapped onto the single im- 
purity Anderson model within DMFT, similar arguments 
must apply. And it is possibly this competition here, that 



FIG. 9: (Color online) The splitting between the up and down 
spin bands is plotted against field for various U values. 
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FIG. 10: (Color online) Comparison of the experimental mag- 
netization in liquid "^He (filled circle) with theory (solid line). 



is ultimately responsible for the first order metamagnetic 
transition. We see that the field induced transition is very 
different from the temperature induced transition, or the 
pure interaction induced Mott transition at zero field. 
In the latter especially, as one approaches Uc2 from the 
metallic side, the Fermi liquid scale Zi* vanishes contin- 
uously, and upon crossing into the Mott insulator side, a 
large and finite gap appears discontinuously in the den- 
sity of states. At Uci the reverse happens, i.e the gap 
in the Mott insulator vanishes continuously on decreas- 
ing [/, and on crossing Uci, a Fermi liquid with a finite 
quasiparticlc weight Z is found. At finite field, however, 
neither the Fermi liquid scale, nor the gap, vanishes at 
the metamagnetic transition. 

Finally we would like to compare our theory with 
experiments. The experimentally measured magnetiza- 
tion^i as a function of effective field at temperatures 
T ^ 70 — 90mK is shown in figure [TUj as filled circles. 
The authors of the experimental paper concluded that 
the first order metamagnetic transition scenario proposed 
by VoUhardti^ cannot be applicable to the case of liquid 
•^Hc, since metamagnetism is absent in the experimental 
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data, although a shght mctamagnctic kink can be distin- 
guished at a field of about 180T. We argue that the appli- 
cability of the half-filled Hubbard model cannot be ruled 
out yet, since as our figure [1] shows, the magnetization 
can be a smooth function of the field for low to inter- 
mediate interaction strengths. According to Vollhardt^, 
U is typically about 15K for ^Hc, while the parameter 
U/Ep ~ 0.8, where Ep is the Fermi energy. So this 
does seem to be an intermediate coupling scenario. In- 
deed, if we superimpose the calculated magnetization for 
U = l.Ot^ onto the experimental data (after a multiplica- 
tive scaling of the field axis), we see surprisingly good 
agreement. The Gutzwiller approach yields a discontinu- 
ous metamagnetic transition for rather weak interactions, 
while we find that metamagnetism does not develop un- 
til a reasonably strong interaction strength U > 1.3t* 
within the IPT+DMFT approach employed here. 



ties. For intermediate U values, we find a first order itin- 
erant metamagnetic transition, in agreement with ear- 
lier theoretical studies. An important benchmark for our 
study was the agreement with the NRG-I-DMFT calcu- 
lations^l regarding the onset of metamagnetism and the 
coexistence region in the h — U plane. We conjecture that 
the first order metamagnetic transition is a result of the 
competition between Kondo screening and Zeeman cou- 
pling. The metamagnetic transition is found to be very 
different in nature from the zero-field interaction driven 
Mott transition. Beyond U = 2.2t^:, we find that even an 
(numerically) infinitesimal field is sufficient to push the 
system from a zero field paramagnetic state to a polarized 
quasi-ferromagnetic state with a substantial magnetiza- 
tion. Finite temperature transport and dynamics studies 
would be instrumental in furthering an understanding of 
this transition. 



IV. CONCLUSIONS 

In this paper, we have studied the effect of magnetic 
field on correlated electron liquids within the framework 
of the half filled Hubbard model using a dynamical mean 
field theory. We compute the the single particle dynam- 
ics, magnetization and other zero temperature proper- 
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